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Editor s Preface 



A generation has passed since the late Emil Artin’s little classic on the 
gamma function appeared in the Hamburger Mathematische Eitizelschriften. 
Since that time, it has b«en read with joy and fascination by many thousands 
of mathematicians and students of mathematics. In the United States (and 
presumably elsewhere as well), it has for many years been hard to find, and 
dog-eared copies and crude photocopies have been passed from hand to hand. 
Professor Artin’s monograph has given many a student his first look at genuine 
analysis — the delicacy of its arguments, the precision of its results. Artin had 
a deep feeling for these aspects of analysis, and he treated them with a master’s 
hand. His undergraduate lectures in the calculus, for example, were filled with 
elegant constructions and theorems which, alas, Artin never had time to put 
into printed form. We may be all the more grateful for this beautiful essay, 
and for its appearance in a new English edition. Various changes made by 
Artin himself have been incorporated in the present edition. In particular a 
small error following formula (59) (this edition) was corrected on the basis of 
a suggestion by Professor Borge Jessen. 

Finally, thanks are due to the translator, Mr. Michael Butler, and to the 
firm of B. G. Teubner for English-language rights. 
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I have written this monograph with the hope of filling in a certain gap 
which has often been felt to exist in the mathematical literature. Despite the 
importance of the gamma function in many different parts of mathematics, 
calculus books often treat this function in a very sketchy and complicated 
fashion. I feel that this monograph will help to show that the gamma function 
can be thought of as one of the elementary functions, and that all of its basic 
properties can be established using elementary methods of the calculus. 

As far as prerequisites are concerned, the reader need only be well acquainted, 
with calculus, including improper integrals. Some of the more important 
concepts needed will even be introduced and discussed again in the first chapter. 
With this background the reader should have no trouble understanding every- 
thing but the later parts of the last two chapters, which do assume some knowl- 
edge of Fourier series. But then, these parts of the monograph can be passed 
over on a first reading without any difficulty whatsoever. 

The following parts of the theory will not be discussed: 

(1) Extension to complex variables. For those familiar with the theory of 
complex variables, it will suffice to point out that for the most part the expressions 
used are analytic, and hence they retain their validity in the complex case 
because of the principle of analytic continuation. The only parts of the theory 
that really need to be changed are those dealing with approximations. This 
certainly should not be much of an obstacle. 

(2) Holder’s theorem showing that the gamma function does not satisfy 
any algebraic differential equation. 

(3) Rummer’s series and the integral representation of log F(x). 

(4) The formula for the logarithmic derivative of F(x). All the necessary 
expressions for this can easily be worked out by the reader. 

I have chosen the integral as my original definition of the gamma function 
because this approach saves us the trouble of proving the convergence of 
Gauss’ product. Any other analytic expression having the characteristic proper- 
ties of the gamma function could just as well have been used. The whole theory 
will then be deduced using the concept of log convexity. This method comes 
from Bohr and Mollerup. 1 

Emil Artin 
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Convex Functions 



Let/(.v) be a real-valued function defined on an open interval a < x < b 
of the real line. For each pair jcj , x 2 of distinct numbers in the interval we form 
the difference quotient 

9>(*i - x i) = " {^ = <p(*2 > * 1 ). ( 1 • 1 ) 

*M x 2 



and for each triple of distinct numbers Xj , x 2 , x 3 the quotient 



uj, x _ • * 3 ) - <fi x 2 . *a) 

, .v 2 , x 3 ) 



(*> - X 2 )f( Xl ) -f (*1 ~ Xs)f(x 2 ) + (X 2 ~ xJfjXz) 
(.V| A'fc) (a*2 *3) (#3 ^1) 



The value of the function ^(x t , x 2 , .v 3 ) does not change when the arguments 
* 1 . *2 » *3 are permuted. 

/(.v) is called convex (on the interval ( a , b)) if, for every number x 3 of our 
interval, <p(x, , .v 3 ) is a monotonically increasing function of Xj . This means, 
of course, that for any pair of numbers .v, > .v 2 distinct from x 3 the inequality 
<p(x x , .v 3 ) > <p(x, , x 3 ) holds; in other words, that '/ / (v 1 , ,v 2 , ,v 3 ) 2? 0. Since the 
value of V is not changed by permuting the arguments, the convexity of f(x) 
is equivalent to the inequality 



n *, ,x 2 ,x 3 )'^o 



( 1 . 3 ) 



for all triples of distinct numbers in our interval. 

Suppose g(x ) is another function that is defined and convex on the same 
interval. By adding (1.3) to the corresponding inequality for^(.v), we can easily 
see that the sum /(.v) + g(x) is also convex. More generally, suppose /,(*), 
/*( x),/ 3 (x) ■■■ is a sequence of functions that are all defined and convex on the 
same interval. Furthermore, suppose that the limit lim„ ^,f n {x) — f(x) exists 
and is finite for all x in the interval. By forming the inequality (1.3) for f„{x) 
with arbitrary but fixed numbers .v, , .v 2 , x 3 , and then taking th£ limit as 
n — ► o°, we see that f(x) is likewise convex. This proves the following theorem: 



I 
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THE GAMMA FUNCTION 



Theorem 1.1 

The sum of convex functions is again convex. The limit function 
of a convergent sequence of convex functions is convex. A convergent 
infinite series whose terms are all convex has a convex sum. 

The last statement of this theorem follows from the fact that each partial 
sum of the series is a convex function and the sum of the series is merely the 
limit of these partial sums. 

We are now going to investigate some important properties of a function 
f(x) defined and convex on the open interval (a, b). For a fixed x 0 in the interval 
let x x range over all numbers > .v 0 and x> range over all numbers < x 0 . We 
have 

9>(*i > *o) > - x oY (>- 4 ) 



If x 2 is kept fixed and x t decreases approaching x 0 , the left side of Eq. (1.4) will 
decrease but always remain greater than the right side. This implies that the 
“right-handed” derivative of f{x) exists; that is to say, the limit 



lim <p(x x , x 0 ) = lim , 

*i>*o ■V-'-'o 



for which we shall use the intuitive notation /'(x 0 + 0). Furthermore, the 
inequality (1.4) also shows that 

f'(x 0 H- 0) y{x, , x 0 ). 

If we let x 2 increase, approaching x 0 , we see that the “left-handed” derivative 
/'(x„ — 0) also exists, and that 

/'(.v 0 + 0) ■>/'(*„ ~ 0). (1.5) 

Given two numbers x 0 < .Vj in our interval, we can choose x 2 , x 3 such that 
x 0 < x 2 < x 3 < .v, . Then 

<K V 2 - v 0 ) < <p( v 3 . x 0 ) = ?>(.v 0 , .v 3 ) < <p(x, , x 3 ) = <p(.v 3 , x,). 

If we let x 2 approach ,v 0 and .v 3 approach .v, , we obtain 

/'(* o + OK/'(.v,-0) for x 0 <x,. (1.6) 



This proves that the one-sided derivatives of a convex function always exist 
and that they satisfy the inequalities (1.5) and (1.6). We shall refer to the pro- 
perties (1.5) and (1.6) by saying that the one-sided derivatives are monotonically 
increasing. 
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In order to show the converse, we must generalize the ordinary mean-value 
theorem to cover the case of functions for which only the one-sided derivatives 
exist. The analogue to Rolle’s theorem is the following: 

Theorem 1.2 

Let /(.v) be a function, defined and continuous on a < x < 6, whose 
one-sided derivatives exist in the open interval a < x < b. Suppose 
fi a ) = /(£)• Then there exists a value £ with a < £ jyb such that one of 
the values /'(£ + 0) and /'(£ — 0) is is 0 and the other ^ 0. 

Proof 

(1) If f(x) takes on its maximum £ in the interior of our interval, then 

h 

is ^ 0 for positive h, ^ 0 for negative h. Taking limits, we get /'(£ -f 0) < 0, 

m- o)>o. 

(2) If the minimum £ is taken on in the interior, we obtain similarly 

/'(* + 0)>0,/'(f-0K0. 

(3) If both maximum and minimum are at a or b, then /(x) is constant, 
f '(x) = 0, and £ can be taken anywhere in the interior. This completes the proof. 

The substitute for the mean-value theorem is the following: 

Theorem 1.3 

Let f(x) be defined and continuous on a ^ x ^ b and have one- 
sided derivatives in the interior. Then there exists a value f in the interior 
such that {/(b) — f(a))j(b — a) lies between /'(f - 0) and /'(f + 0). 
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THE GAMMA FUNCTION 



We are now in a position to prove the desired converse. Let f(x) be a func- 
tion defined on the open interval a < x < b. Suppose /(x) has one-sided deriva- 
tives that are monotonically increasing. We contend that f(x) is convex. 

Let .Vj , x 2 , x 3 be distinct numbers in our interval. Since the value of V 7 
does not change under permutation of the subscripts, we may assume that 
.Vj < .v 3 < x, . According to the mean-value theorem, we Can find £, rj with 
.v, <rj <x 3 < £ <x x such that <p(x, , x 3 ) lies between /'{£ — 0) and f(£ + 0), 
and <p(x 2 , x 3 ) between /'(rj — 0) and /'(>? + 0). Therefore (1.5) implies that 

?(*i - *3) >f\t ~ 0) and <p(* s , * 3 ) < f'(r, + 0). 



From Eq. (1.2) we obtain 

'f , (x l . x, , x 3 ) 



f\i - Q) -f(v + 0) 
*1 - *2 



Finally we conclude from (1.6) that 



^(a'i , x 2 , x 3 ) 0, 

and this is the contention. 



Theorem 1.4 

/(.y) is an convex function if, and only if,/(.xr) has monotonically increas- 
ing one-sided derivatives. 

Corollary 

Let f(x) be a twice differentiable function. Then /(x) is convex if, 
and only if, f"{x) 0 for all x of our interval. 

Proof 

/'(x) is monotonically increasing if, and only if, /"(x) 0. 

We now return to Eq. (1.2) and select for x 3 the midpoint (x, -4- x 2 )/2 of 
.v, and x 2 . Assuming for a moment that ,v 2 < ,v l , we have 

■ v 3 x t = V 1 V 3 — 2 (* v i ~ *2)- 

The numerator of x P(x l , .v 2 , .y 3 ) becomes 

(*i - *2) (i/(*i) -r £/(**) -/(x 3 )), 

and the denominator is positive. For a convex function we obtain the inequality 

/ (— L ^ L: —) C i(/(*i) + /(*a)), 



(1.7) 
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which is symmetric in Xj and x 2 and therefore also holds for Xj < x 2 . For 
x t = x 2 it is t rival. 

We shall call a function defined on an interval weakly convex if it satisfies 
the inequality (1.7) for all x l , x 2 of the interval. It is obvious that the sum of 
two weakly convex functions, both defined on the same interval, is again weakly 
convex. It is also obvious that the limit function of a sequence of weakly convex 
functions, all defined on the same interval, is weakly convex. 

Let /(x) be weakly convex. The inequality (1.7) can be generalized to 

/(“ ") ^ ~(/(*i) + /(**) + ••• +/(*„)). (1.8) 



Proof 

(1) We first show that if (1.8) holds for a certain integer n, then it also 
holds for In. Indeed, suppose Xj , x 2 , •••, x 2n are numbers in our interval. 
Replacing x, and x 2 in Eq. (1.7) by 

*1 + ••• + an( j *n+l + ••• + 

n n 

respectively, we have 

/( * 1+ 2 „ +Xz ’‘ ) < i (/(“——) +/{ Xn+1 + „"~ —)) • 

Applying the inequality (1.8) to both terms on the right-hand side, we get the 
desired formula 

/( — " 2n +Xin ) ^Tn U{X ' ] +f{Xi) + "• +/(**«»• 

(2) Next w'e show that if (1.8) holds for n + 1, then it also holds for n. 

W'ith n numbers (x t , x 2 , x„) the number 

= ~ (*1 + ' " + *n) 

also belongs to our interval. If (1.8) holds for n + 1, then 



t nx n+- 1 + *n+l\ 




+ + A *n+l\ 


\ n+ 1 > 




n + 1 ) 



< ~7j-(/(*i) + ••• +/(*«) -r /(* n+l))* 



Transposing the term 1 /(« + l)/(x n+1 ) to the left side, w>e obtain (1.8) for the n 
given numbers. 
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THE GAMMA FUNCTION 



(3) We now combine steps (I) and (2) to attain the desired result. If 
(1.8) holds for any integer n, then step (2) implies that it also holds for all smaller 
integers. Because of step (1) the contention is true for arbitrarily large integers. 
Therefore it must be true for all n. This completes the proof. 

We wish to prove the following theorem: 

Theorem 1.5 

A function is convex if, and only if, it is continuous and weakly con- 
vex. 



Proof 

(1) A convex function is continuous since it has one-sided derivatives. 
It is also weakly convex, as has already been shown. 

(2) Suppose that f(x) is weakly convex, that there are x 2 < x l numbers 
in our interval, and that 0 < p < n are two arbitrary integers. Apply (1.8) 
to the case where p of the n numbers have the value x, and the remaining 
n — p numbers have the value x 2 ■ We obtain 

/ (ir *« + ( 1 1 - t < £/(*») + ( 1 _ t) /{Xi) - {1 - 9) 

Assume now that f(x) is continuous and let t be any real number such that 
0 < / < 1 . Select a sequence of rational numbers between 0 and 1 that converges 
to t. Every term of this sequence is of the form pjn for suitable integers p and n; 
therefore Eq. (1.9) can be applied. Since f(x) is continuous, we can go to the 
limit. We obtain 

/(**, i + ( 1 - t) *a) ^ </(*.) + ( 1 ~ 0/(*i>- (110) 

For any distinct numbers (x\ , x 2 , x 3 ) in our interval we must show that 
0(xj , jt 2 , x 3 ) 0. Since 4> is symmetric, we may assume that x 2 < x 3 < .Vj . 

The denominator of Eq. (1.2) is positive. 

We set t — (,v 3 — x i )j(x l — ,v 2 ); then 

0 < f < 1, 1 — / = — * 

X\ ~ - v 2 

and 

tX, + (1 -t)x 2 = ( x 3 - x 2 ) x v -t- (x ir -x 3 ) x 2 = ^ 

X ^ ~ Xf^ 



Hence Eq. (1.10) implies that 



f( x i) ^ 






which shows that the numerator of <p i9 ^ 0. This completes the proof. 
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Numerous inequalities useful in analysis can be obtained from Eq. (1.8) 
by a suitable choice for /(*). As an example, consider f(x) = — log x for x > 0. 
We have f'\x) = I/x 2 and our function is convex. Therefore Eq. (1.8) implies 
that 

- log ( + n ' ~ l ~ (*og x 3 -|- log x 2 + ••• + log *„), 

hence 

log > log <y x l + ••• +x:„ , 




We now introduce an important concept closely related to that of convexity. 
A function f(x) defined and positive on a certain interval is called log convex 
(weakly log convex) if the function log f(x) is convex (weakly convex). The con- 
dition that f(x) be positive is obviously necessary, for otherwise the function 
log f(x) could not be formed. As an immediate consequence of our previous 
results, we have the following: 



Theorem 1.6 

A product of log-convex (weakly log-convex) functions is again 
log convex (weakly log convex). A convergent sequence of log-convex 
weakly log-convex) functions has a log-convex (weakly log-convcx) limit 
function, provided the limit is positive. 

Instead of the condition that the limit function be positive, we could 
require that the sequence of the logarithms of the individual terms be con- 
vergent. 



Theorem 1.7 

Suppose f(x) is a twice differentiable function. If the inequalities 

/(*)> 0, /(*)/"(*) -(/W>0 

hold, then f(x) is log convex. 

This theorem follows directly from the fact that the second derivative of 
log f(x) has the value 

/(*)/"(*) - (/ X *)) 2 
(/(*))* 
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THE GAMMA FUNCTION 



The properties of log-convex functions mentioned thus far are all more or 
less immediate consequences of the definition. The following remarkable 
theorem, however, is a much deeper result: 

Theorem 1.8 

Suppose /(*) and g(x) are functions, defined on a common interval. 

If both are weakly log convex, then their sum /(*) + g( x ) is also weakly 
log convex. If both are log convex, then f(x) + g(x) is log convex. 

Proof 

It suffices to prove the first statement. The second then follows immediately 
with the addition of continuity. 

Both/(x) and g(x) are positive. If x t , x 2 are numbers in our interval, then 

(/( 2 )) 2 ) and (? (— 2 ** )) <g( x i)8( x t)- 

We have to show that 

(/ (— + g f— < (M) + <?(*,)) (/(*,) g(x 2 )). 

In other words the proof of our theorem amounts to showing that if a x , b x , c x , 
« 2 . , c 2 are positive real numbers with a x c x 0 and a 2 c 2 — b\ ^ 0, 

then 

(«i + <* 2 ) (<h + c 2 ) - ( b x + b 2 f > 0. 

Consider the quadratic form a,.v 2 + 2 b x xy + c x y 2 where a x > 0. We have 

-f- 2 b t xy + c x y 2 ) =-. (a x .x + b x yf + (a x c x - b\) y 2 . 

If a x c x — b\ > 0, the quadratic form never takes on a negative value, what- 
ever x, y may be. On the other hand, if a x c x - b' 2 < 0, the quadratic form takes 
on the negative value a x c x - - 6 2 for y = 1 , x — - (b x !a x ). 

Our conditions imply that neither 

a x x 2 -r 2 b x xy ~f c x y 2 nor a^x 2 + 2b.yxy -j- c 2 y 2 
takes on negative values. Therefore 

( a i + “ 2 ) A' 2 + 2 (b x -f- b 2 ) xy + (cj + c 2 )y 2 

will not take on negative values. Consequently 



(«i + a 2 ) (c, + c 2 ) - (b x + b 2 f 0. 
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This completes the proof of the theorem. 

The reader who enjoys working with identities can check the validity of 
the following, for an alternate proof: 

<WK + (O + O) - (*i + KY) 

= a 2 (a x + a 2 ) ( a x c x - b\) + a x (a x + a 2 ) (a 2 c 2 - 6 2 ) + (a x b 2 - aj) x f. 

If 

U] ^ 0, o 2 ^ 0, ^ 0, and a 2 c 2 b 2 ^ 0, 

the right side is ^ 0 and the conclusion follows. 

Other important facts can be obtained by combining our previous results. 
Suppose f(t, x) is a function of the two variables x and t, which is defined and 
continuous for t in the interval a < t < b and x in some arbitrary interval. 
Furthermore, for any fixed value of t, suppose that /(t, x) is a log-convex, twice 
differentiable function of x. For every integer n we can build the function 

F n (x) — h(f(a, x) +/(a + h, x) + f (a + 2h, x) + + f(a + (n - l)h,x)), 

where h = (b — a)/». Being the sum of log-convex functions, F„(x) is also 
log convex. As n approaches infinity, the functions F„(x) converge to the integral 

f f(t, x) dt; 

J a 

hence this integral is also log convex. 

Suppose that f(t,x) only satisfies our conditions in the interior of the 
t interval, or that the upper bound of the interval is infinite. If the improper 
integral 

f /(<>*)<* 

J a 

exists, then it is log convex. This follows directly from the fact that an improper 
integral is the limit of proper integrals over subintervals. Hence, as the limit 
function of log-convex functions, it is also log convex. 

In this book we will only have to test integrals of the form 

r b 

<p(t) t x 1 dt 

** O 

for log convexity, where <p(t) is a positive continuous function in the interior 
of the integration interval. If we take the logarithm of the integrand and then 
differentiate twice with respect to x, we get 0. 
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Theorem 1.9 

If 95(f) is a positive continuous function defined on the interior of 
the integration interval, then 

|%(<) t x ~ l dt 

J (t 

is a log-convex function of x for every interval on which the proper or 
improper integral exists. 

The following theorem is quite obvious: 

Theorem 1.10 

If f(x) is log convex on a certain interval, and if c is any real number 
# 0, then both the functions f(x + c) and f(cx) are log convex on the corres- 
ponding intervals. 



[ 2 ] 



The Euler Integrals and the 
Gauss Product Formula 



The theory of the gamma function was developed in connection with the 
problem of generalizing the factorial function of the natural numbers, that is, 
the problem of finding an expression that has the value n\ for positive integers n, 
and that can be extended to arbitrary real numbers at the same time. In looking 
for such an expression, we come upon the following well-known improper 
integral: 

j e~‘ t n dt = n\ 

J 0 

This suggests replacing the integer n on the left side by an arbitrary real number 
(provided the integral still converges) and defining x! for arbitrary x as the value 
of this integral. Rather than doing precisely that, we will follow the custom 
of introducing a function that has the value (w — 1)! for positive integers n. 
Namely, 



,CC 

r( X ) = I 

■> 0 



We still must determine the values of x for which this integral converges. The 
integrand is smaller than t x ~ l when t is positive; therefore 



I e~‘ t dt < ( t *- 

J < • < 



1 € X 

1 dt= , 



For x > 0, 



| e~* t < 



is bounded from above by 1/x. If we hold x fixed and let e decrease, the value 
of the integral increase monotonically. This means that 



f e~‘ t x ~ l dt = lim f e~‘ t x ~ x dt 
Jo «-o 3. 



exists for all positive x. 
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THE GAMMA FUNCTION 



When t is positive, every term of the series for e‘ is positive, and 
the inequality e ( > t n jn\ holds for all integers n. Hence e~ l < n\jt n , which gives 
another inequality for the integrand, namely, e- l t x ~ l < n!/t n+1_a: . Therefore, by 
holding x fixed and choosing n > * + I, we can make «!/(« — a) an upper 
bound for 



ft 

e-U*- 1 , 
J 1 



But the value of this integral increases as 8 increases, and thus 



,» ..a 

e~‘ t x ~ l dt = lim I e~ l t x ~ l dt 
J , a -mc J t 



exists. This implies that our definition, Eq. (2.1), is meaningful for all positive 
real x. 

If we replace x by x + 1 in Eq. (2.1) and integrate the approximating 
integral by parts, we get 

r a a ,a 

e~ l t x dt = - e~* t x + * e~* t x ~ 1 dt 

» < » « 

.■a 

= e~'t z — e~*8 x + x e~ l t 1 " 1 dt. 

* « 

The term is again smaller than n\j8 n ~ x . As e approaches 0 and 8 approaches 
+ <», the first two terms on the right side vanish, and we have the formula 

f(x + 1) = x r(x). (2.2) 



This functional equation is basic for the development of the rest of the theory. 
It represents a generalization of the identity n\ = n(n — 1)! for nonintegral 
values of n. Suppose the value of the gamma function is known on the interval 
0 < x ^ 1 . With the help of Eq. (2.2) we can easily calculate its value on the 
interval 1 < x < 2, then again on the next interval of length 1, and so on. By 
repeated application of Eq. (2.2), we get 

f(x + n) = (x + n — 1) (x + n — 2) ••• {x + 1) xf(x) (2.3) 



for every positive integer n. 

Equation (2.1) is only a definition for positive x. Now we want to extend 
this definition to include negative real numbers. If x lies in the interval 
— n < x < — n+1, we define the value of the gamma function at x by 
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If x is a negative integer orO, the right side of Eq. (2.4) is meaningless. We will 
consider F(x) as undefined for these particular numbers. Otherwise, the left 
side of Eq. (2.4) is well defined, since the argument (x + n) on the right 
lies in the interval 0 to 1 . This extended definition is obviously so constructed 
that the functional equation, Eq. (2.2), always holds. 

The two paragraphs above clearly show that Eq. (2.2) does not determine 
the gamma function uniquely. If f(x) is any arbitrary function defined on the 
interval 0 < x ^ 1 , we can set 

/(* + ») = (* + «— 1) (*+ n - 2) (x + 1) */(*) (2.5) 



for 0 < 1, and 



/(*) = 



x(x + 1) ••• (* + n 



-J(x + ») 



for — n < x < — n + 1. Thus f(x) is defined for all real numbers, with the 
exception of 0 and the negative integers, in such a way that the functional 
equation f(x + 1) = xf(x) always holds. This certainly makes our definitions, 
Eqs. (2. 1 ) and (2.4), seem rather arbitrary. If we keep our original problem in 
mind, it is quite natural to want Eq. (2.2) to hold. It is the appropriate generaliza- 
tion of an elementary property of the factorial function. But an infinite number of 
arbitrary functions can be found that share this property with the gamma 
function. What singles out 7'(x) from all the other possible functions we could 
have defined ? One glance at the integral in Eq. (2.1) shows that 



n i) = i. 



and therefore 

m = («-!)!• ( 2 . 6 ) 



Furthermore, P(x) is continuous and differentiable. (This will be proved later on.) 
But even so, an infinite number of other functions that also have these properties 
can be found. 

Our integral in Eq. (2.1), however, has another property that catches the eye. 
It is log convex. This fact follows immediately from our conclusions concerning 
integrals and log convexity in Chapter 1. Intuitively, it means that the curve 
y = log r(x) is very smooth. Strange as it may seem, this property is enough 
to single out T(x) from all the other solutions of the functional equation 
f(x + 1) = */(*). 

We shall now prove the following: 
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Theorem 2.1 

If a function f(x) satisfies the following three conditions, then it is 
identical in its domain of definition with the gamma function: 

(1) /(*+ l)=x/(x). 

(2) The domain of definition of f(x) contains all x > 0, and is log 
convex for these x. 

(3) /(1) = 1. 



The existence of a function with these properties (the gamma function) 
has already been proved. 

Suppose /(x) is a function that satisfies our three conditions. Then Eq. (2.5) 
is valid because of condition (1), and /(«) =(« — 1)! for all integers n >0 
because of condition (3). It suffices to show that f(x) agrees with T(x) on the 
interval 0 < x ^ 1 . If this is the case, then /( x) must agree with .T(x) every- 
where because of condition (1). Let t be a real number, 0 < x ^ 1, and n 
an integer ^ 2. The inequality 

log/(- 1 ~n) - log f(n) log /(x + n) - log/(w) log/(l + n) - log/(n) 



( - 1 + n) - n 



(* + n) — n 



(1 + n) - n 



expresses the monotonic growth of the difference quotient for particular values, 
and is therefore valid because of condition (2). Since /(«) = (n — 1)1, we have 

i , , . . log/(x + n) — log (n — 1)! _ , 

log (n - 1)< — ^ — - — ' ^ log n 



log ( n — 1)* ( n — 1)! < log f(x + w) < log n x (n — 1)!. 

But the logarithm is a monotonic function; hence 

(n - 1 ) x (n - 1)! </(x + n) sg n x {n - 1)!. 

With the help of Eq. (2.5), we get the following inequality for f(x) itself: 

(n - !)*(« - 1)1 n x {n - 1)! 

x(x + 1) ••• (x + n - 1) ^ x{x + 1) ••• (x + n - 1) 

_ n x n\ x + n 

x(x + 1) ••• (x + n) n 

Since this inequality holds for all n ^ 2, we can replace n by n + 1 on the left 
side. Thus 

«^*1 < ft \ * + » 

x(x + 1) (x + n) ^ x(jc + 1) ••• (ac + n) n 
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An easy calculation gives the inequality 

* x + n ^ x(x + 1) •••(* + ») "" ^ X " 

As n approaches infinity, we get 

n l i2 x(x + 1) •" {x + «) ‘ 

But F(x) is also a function that satisfies our three conditions. Hence the relation 
we have just derived is still valid if we put r(x) instead of f(x) on the left side. 
This completes the proof of the theorem. 



As a corollary we have the formula 



Hx) = lim t • 

V ’ n->a: x(x + 1) •" (X + «) 



Actually Eq. (2.7) was only proved for the interval 0 < x < 1. To show that it 
holds in general, we denote the function under the limit sign by r n (x). It is 
easy to see that 

r n {x + i) = x r n (x) x + l +l , W = y x + n n +l r n (x + i). 



These two expressions help clarify the following fact: As n approaches infinity, 
if the limit in Eq. (2.7) exists for a number x, it also exists for x + 1. Conver- 
sely, if it exists for x + 1 and x # 0, it also exists for x. Hence the limit 
exists for exactly those values of x for which T(x) is defined. If we denote the 
limit in Eq. (2.7) by /(x), we get the equation /(x + 1) = x/(x). Since /(x) 
already agrees with T(x) on the interval 0 < x < 1, it must also agree every- 
where else. Equation (2.7) was derived by Gauss, and it is often used as the 
fundamental definition of the gamma function. 

Another form of Eq. (2.7) which is important in the theory of functions 
was derived by Weierstrass. A simple manipulation shows that 



r n {x) 




g x! 1 gxli 

1 + x/1 1 + x/2 



e xln 

1 + Xfn ■ 




But the limit 
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exists.* It is often called Euler’s constant. Therefore, we have 



1 e xli 1 ^ 

r(x) = e-c* - lim n TTT = «" c * - EL - (2-8) 

X *•-»« “ 1 4- xjt X 1 + x/l ' 



where II is the product symbol. 

We shall now show that the function /’(x) can be differentiated as often as we 
please; that is, T(x) has derivatives of arbitrarily high order. Because of the 
functional equation (2.2), it suffices to prove the assertion for positive x. But 
r(x) > 0 for x > 0, therefore log T(x) is defined. From (2.8) we get 



log r(x) = - Cx - log x + lim £(*_ log (l + y)) 



= - Cx - log x + j? ( 4 - - log (l + 4 -)) . (2.9) 

We will now proceed to prove our assertion for the function log T(x). The 
conclusion for T(x) will then follow immediately from r(x) = s l0 * r(l1 . Can we 
take the derivative of the series in Eq. (2.9) term by term ? This can be done 
if the new series — the one obtained by termwise differentiation — is uniformly 
convergent. Should this be the case, the left side of Eq. (2.9) is differentiable. 
The differentiated series obtained from Eq. (2.9) is 



* To prove this we can set 

C„ = ^ *°K" and D« = C„ 

1 2 n 

which gives us 

Cn+i — C u = — — — log (l + — | , 
n + 1 V n l 

and 



£>.♦1 — £>„ = - tog (l + — ) . 

n \ n / 



The elementary inequality 



I 

n + 1 



(■ 4 t) 




n 



shows that the sequence C„ decreases monotonically, whereas the D n increases monotonic* 
ally. Furthermore 2), < C,; hence, £>, = 0 is a lower bound for the C„. In other words, 
the sequence C„ converges to a limit. 
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Because x is positive, the general term of this series is smaller than x/i*. If we 
restrict x to an arbitrary interval 0 < x < r, then the general term is smaller 
than r/ «*. This number is completely independent of x. The series 



« r /111 \ 

t 1 * T ( 1* + 2 2 + 3* + " 



converges; therefore, our differentiated series converges uniformly. But this 
means log T(x) is differentiable, and that, on the interval in question, 



d r\x) 

Tx 108 r( ') = Tw 



=-c-| + £(i- r y. 



But the choice of this interval was arbitrary, which implies that Eq. (2.10) 
holds for all x > 0. 

Suppose we take the derivative once again. We get the series 



± + v_ !— . 

X* ^ H (x + o* 



** - + 

The general term is smaller than l/i* because x > 0. This series obviously 
converges uniformly for positive x. Repeated differentiation leads to ever better 
converging series, all of which converge uniformly for positive x. This shows 
that log jT(x) can be differentiated as often as we please. We have the formula 

d*-'/r(x)\ A(-1)*(A-1) 

(* + »•)* ’ k > 2 - ( 2 ’ n > 

There is no trouble in extending the validity of Eq. (2.1 1) to include negative x. 
This is easily done with the help of Eq. (2.2). We merely determine the functional 
equation that the left side of Eq. (2.11) satisfies, and then show that the right 
side also satisfies it. Both these steps are quite obvious. 

The case k — 2 is of special interest. The function 

d ( r\x) \ 



L ( \ 

dx \ F(x) ) 



is always positive; therefore, the following inequality holds for all x: 

T(x) r\x) - (T'(x)) 2 >0 or T(x) T"(x) > (/»)* ^ 0. 

This shows that the functions T(x) and -T"(x) are either both positive or both 
negative for each particular value of x. Consequently the function | f(x) | 
is convex. When is /’(x) positive and when is it negative ? We already know that 
T(x) is positive for positive x. It follows from Eq. (2.4) that /’(x) has the sign 
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(— 1) K on the interval — m < x < — n + 1. Furthermore, Eq. (2.4) shows 
that r(x) has a very large absolute value whenever x gets close to zero or a 
negative integer. If the well-known values of r(x) for the positive integers are 
also taken into account, we get a good idea of what the curve y = l n (.v) looks 
like. (The reader is encouraged to draw a sketch.) 

The integral of Eq. (2.1) is due to Euler, and it is referred to as Euler's 
second integral. He also discovered another integral related to the gamma 
function, which is called Euler's first integral: 

B( x, y) = f 1 /*->( 1 - t)*- 1 dt. (2. 1 2) 

J o 

This time we have a function of the two variables x and y. We want to prove 
that the integral exists whenever x and y are positive. First, we write our integral 
as the sum of two integrals, one from 0 to and the other from to 1. The 
integrand of the first is always smaller than t I_1 (l — f)- 1 and hence smaller 
than 2 t x_1 . In the second the integrand is always smaller than < _1 (1 — t) v_1 
and hence smaller than 2(1 — t) v ~ x . The method we used to prove the existence 
of Euler’s second integral can now be applied to these two integrals. The rest 
of the details are left to the reader. 

If we replace x by x + 1 in Eq. (2.12) and write the integral in the form 

B(x + \,y) = J‘(i - o x+v -* dt < 

we can integrate by parts. We get 

• 1-d / t vX 

J, c-r-H* 

il — tY +i/ , t \ x i l_a r l ~ d x it \ 1-1 I 

■ - “Try - (rb) I. + /. rb<‘ - ,) **’ (rh) (rb? * 

_ (!_-.) ». •->'(! -8) ; f'-* _ , r , 

x+y x+y) f 



(1 - t) v t* — 8 V (1 -8) 

x -t- y 



If we let e and 8 converge to zero, we get the following functional equation: 

«<* + !->-) = y)- 

Now we hold y fixed and consider the integral of Eq. (2.12) as a function of x. 
In order to obtain a function that satisfies the functional equation (2.2) we set 

/(*) = B(x, y) T(x + y). 

This function obviously satisfies condition (1) in Theorem 2.1. Furthermore, 
/( x) is the product of two log-convex functions and is therefore log convex 
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itself. The log convexity of £(x,_y), regarded as a function of x, follows imme- 
diately from our theorem on the log convexity of integrals. T(x + y) is obviously 
log convex. This means that /(x) also satisfies condition (2) in Theorem 2.1. 
Condition (3), however, does not hold. We have 

£(!.>’)= |‘ (l - O'" 1 *=-7’ 

j o y 

and therefore 

/(l) = lr(l +y) = r(y). 

But this is not really a serious difficulty. Given any function £(x) that satisfies 
conditions (1) and (2), we can always construct a function that also satisfies (3). 
Condition (2) implies that £(1) >0, which means we can form the quotient 
g(x)lg(l). This function satisfies all three conditions and therefore is E(x). In 
other words, we have 

S( x ) = g{ 1) A*)- 

In our particular case we get 

/(x) = r(y) r( X ). 

But this means we have evaluated the integral in Eq. (2.12): 

(2.13) 

This formula holds for all positive x and y. 

By setting x = £ and y — £ in Eq. (2.13), we get an integral that can 
easily be evaluated. The substitution t — sin 2 <p gives us 

m)y = 2f rl d<p = n. 

* 0 

But rjt is positive; therefore, we have the following remarkable and important 
identity: 

r(i) = . (2.14) 

Using Eq. (2.14) and the functional equation (2.1), we can easily calculate the 
value of P(n + £) for integral n. 



[ 3 ] 



Large Values of x and the 
Multiplication Formula 



Can we find an elementary function that gives an accurate approximation 
of r(x) for large values of x ? If the growth of ft! is estimated, it is found to 
increase with n faster than rt n e~ n , but not quite as fast as n n+1 e~ n * In other 
words, the growth of T(«) is caught between n”~ l e~ n and n n e~ n . This suggests 
that we consider a function of the form 

/(*) = **-»' 2 (3.1) 

in order to study the behavior of f(*) for large x. Our goal is to make /(*) satisfy 
the basic conditions for the gamma function by choosing fx(x) in an appropriate 

way. 

If we replace x by x + 1 in Eq. (3.1) and divide the resulting expression 
by Eq. (3.1), we get 

f( X *) _ /] i J\ X+1/2 ^(x+U-.uU) 

/(*) V x! 

This shows that f(x) satisfies condition (1) in Theorem 2.1 if, and only if, 

p(x) - h (x +!)-(* + *) log (l + !■) - 1, (3.2) 

holds for n(x). 

* If we consider the elementary inequalities 

for k = 1,2, ••• , (»i — 1), and multiply them together, we get 

n*' 1 , n" 

— < e "' 1 < . 

(ft— 1)! (ft— D! 

This leads to the approximation 

Cfi"e - " < n! < en''* , e~". 
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We denote the right side of Eq. (3.2) by j-(*). A function n(x) with this 
property is easy to find. If we set 

00 

m(*) = X (3.3) 

«-o 



then Eq. (3.2) holds, provided the infinite series in Eq. (3.3) converges. Let us 
postpone the proof of convergence for a moment and consider condition (2) 
of theorem 2.1. 

The factor x x ~ ln e~ x in Eq. (3.1 ) is log convex because the second derivative 
of its logarithm, Ijx + | * 2 , is always positive when x is positive. If we can 
show that the factor e> x{x) is log convex, in other words that n(x) is convex, then 
f(x) also satisfies condition (2). This means that the function f(x) determined 
by the particular n(x) we defined in Eq. (3.3) will agree with F(x) to within a 
constant factor. Our n(x) is convex if the general term of the series g(x + n) 
is convex. To show this, it suffices to prove the convexity of g(x) itself. But 
we have 



*"(*) 



1 

2* 2 (* + l) 2 



> 0 . 



The convergence of the series in Eq. (3.3) still remains to be shown. We 
will combine this with an approximation of the function /*(*). Let us begin by 
considering the expansion 






which is valid for \y | < 1. Now we replace y by 1/(2* + 1). The resulting 
expansion is valid for positive * because 1/(2* + 1) < 1 whenever * > 0. 
We multiply this equation by 2* + 1 and bring the first term on the right side 
over to the left side: 



(* + £)log(l +1) - 1 =g(x) 

1 1 1 

3(2* + l) 2 1 5(2* + \)* + 7(2* + 1)» + ’ 

This expression again shows that g(x) is convex, since every term on the right 
side is convex. Now we can approximate £(*). If the integers 5, 7, 9, ••• are all 
replaced by 3, then the value of the right side increases. The result is an infinite 
geometric series, having l/(3(2* + l) 2 ) as its first term and 1/(2* + l) 2 as its 
ratio. Its sum is 



II II 

3(2* + l) 2 1 - (1/(2* + l) 2 ) 12*(* + 1) “ 12* 



1 

12(* + 1) • 
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But g(x) is positive, hence 

0 < g( x ) < Y2* — 12(.v + 1) • 

Since every term of the series in Eq. (3.3) is positive, it suffices to show the 
convergence of 






"U2f* + ») \2(x + n + 1) / ’ 

which converges trivially to the limit ]/l2x. This not only proves our assertion, 
it also gives the approximation 



In other words, 



0 < n(x) < 



*<*> - js 



where d is a number independent of x between 0 and 1. 

By a suitable choice of the constant a, we get 

r(x) = ax x - 1 ' 2 e - x+ f lx> = ax x ~ 112 e ~ x+en2x . (3.4) 

If we let x be an integer n and multiply the expression by n, we get the appro- 
ximation 

«! = an n+l/2 e~ n+6,l2n . (3.5) 

We are now going to find the exact value of this constant a and determine 
some other important constants at the same time. 

Let p be a positive integer. We consider the function 

/w _,.r(£)r(iii)...r(i ± ^>), 

for x > 0. The second derivative of log p x is zero, and each of the functions 
r((x + i)jp) is obviously log convex. This implies that f(x) is also log convex. 
If we replace x by x + 1 , p x takes on the factor p, r((x + i)!p) goes over into 
the next factor, and r({x + p — 1 )jp) becomes 

r (i + ') = 7 r (i)- 

In other words, f(x) is multiplied by x. Our function again satisfies the conditions 
(1) and (2) in Theorem 2.1; therefore, 
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where a v is a constant depending on p. For x = 1 in Eq. (3.6), we have 

».-f r (!) r (7)"' r (£)- < 3 - 7 > 

If we set x = k\p in Eq. (2.7), then a simple manipulation gives 

r (A) _ r n l "‘n\p n ^ 

l p ) k(k +p)(k + 2p)---(k + np) • 



Now we set k = 1, 2, •••,/>, one after the other, and multiply all these expressions 
together. Factors of the form ( k + hp ) appear in the denominator, where k 
runs from 1 to p, and h runs from 0 to n. For h = 0 we get the numbers from 
I to p\ for h — 1 , the numbers from p + 1 to 2 p; and so on. The product in 
the denominator is obviously (np + p)\. The final result is 



a„ = p lim 



n u,+l) I2^ n \y pnv+\ 



n-xo (np -f p)\ 



The well-known infinite product 



1 — lim ( 1 + A) (l +—*)'•• (l + —) > 

oo \ tip! \ np) \ tip) 



which can be written as 



l — lim ( np+ ?l 1, 

«-« (npy. (npy ’ 



can now be applied. If we multiply this last expression with the above identity 
for a v , we obtain 

.. (n!)" p nv 

a„ = p lim t . 
p r n^co (np)\ 

But Eq. (3.5) implies that 

(k!) p = a v ri' v *“ n e~ nv e 9lVinn , 

(np)\ = a(np) nv+112 e~" p e ^ lnnv . 

After making the appropriate substitutions above, we obtain 



a p ~ y/P aP l lim e 



(0 1 p/12n)-(0 ! /12ni.) 



and finally 



a v — Vp a 1 ’-' 1 . 
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By evaluating a t with the help of Eq. (3.7) and then comparing the result with 
Eq. (3.8), we get 

<2 2 =2r(£)r(l) = 2\/*=«V2. 

But this determines the exact values of our constants: 

a = y/l-n and a P — /> 1/2 (2rr) ,p ~ 1,/s . 

Now we gather together all the important expressions from this, chapter: 



r{x) = V2n x 1 - 1 ' 8 e-* +p < x ', 

M*) = X (* + n + i) lo § ( ! + - 1 = ^ - 0 < 0 < 1. 



w! = V2»r h" + i/2 e~ n+t/itn . 



(3.9) 






In particular, for p = 2 






The formulas in Eq. (3.9), which describe the behavior of F(x) for large 
values of x, are called Stirling's formulas. If our approximation of n(x) is used, 
the accuracy of the formula for r(x) will increase as * increases. This is also 
true for estimates of n! The relative accuracy for n ^ 10 is already quite high. 

The functional equation (3.10), discovered by Gauss, is called Gauss’ 
multiplication formula. By replacing x by px in Eq. (3.10), we obtain an expres- 
sion for r(px) as the product of factors, each of the form /’(* + (kip)). This fact 
gave rise to the name “multiplication formula.” The most important special 
case is p = 2. It was discovered by Legendre and is often referred to as Legendre’s 
relation. 



I 4 ] 

The Connection with sin x 



The gamma function satisfies another very important functional equation. 
In order to derive it, we set 

9 ?(x) = r(x) r( 1 — x) sin nx. (4.1) 

This function is only defined for nonintegral arguments. If we replace x by 
x + 1, then r(x) becomes xT(x). The function T(I — x) becomes 

and sin nx changes its sign. This means that <p(*) is left fixed, and is therefore 
periodic of period I : 

9>(* + 1) = ?(*)• ( 4 - 2 ) 

The Legendre relation can be written in the form 

r(y)r(-tl) =w-r ( ^). 

where b is a constant. Actually, the exact value of b was determined in Chapter 3. 
But this extra information need not (and will not) be assumed here. As far as 
we are concerned now, b is just some particular constant. 

In the expression above, we replace x by 1 — x: 

r (4r) r (> - t) - b2x ~ l r v ~ x) - 



Now we consider 



* (y) - f-T 1 ) - r (f > ' p (> - t) si " f r (H- 1 ) r ( L ? ! ) - T 



= — P(x) r(l — x) sin TTXy 



and we get the relation 
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where d is a constant depending on b. The exact value of d is not important 
here. 

Since both f , (.v) and sin x can be differentiated as often as we please, <p(x) 
also has this property. Because of the functional equation (2.2), we can write 

, , r(l + x) . 

<p( x) = — - r(l — x) sin nx 

— T(1 + a) .T(l — x) ^ — | — — . 

where the power series converges for all values of x. But the right side of this 
equation is also defined for x = 0, and it represents a function having derivatives 
of all orders at this point. This suggests that we extend our definition, Eq. (4.1), 
by giving <p(x) the value n at x — 0. Because our function is periodic, we define 
n to be the value of <p(x) for all integral arguments. <p(x) is now continuous 
everywhere and has derivatives of all orders at every point. The relation of 
Eq. (4.3) was only proven for nonintegral x. But if we let x approach an arbitrary 
integer, the validity of Eq. (4.3) for all a- follows from continuity. As far as the 
sign of <p(x) is concerned, Eq. (4.1) shows that <p(x) is positive on the interval 
0 < x < 1. Because of Eq. (4.2), this is now true for all x. 

Our goal is to prove that y>(x) is a constant. Let g(x) denote the second 
derivative of log<p(x). g(x) is also periodic of period 1. Because of Eq. (4.3), it 
satisfies the functional equation 

1 (f (y) + g (— y-)) = g(x). (4.4) 

Since g(x) is continuous on the interval 0 ^ x < 1 , it is bounded on this interval, 
sa y> I g{ x ) I ^ M. But this inequality holds for all x, because g(x) is periodic. 

The following argument shows that g(x) vanishes. Equation (4.4) gives us 
the inequality 




This means that the upper bound can be pushed down from M to M/2. If we 
repeat the process again, we get M/4 as an upper bound, and so on. In other 
words, the upper bound for g(x) can be made as small as we please. This implies 
that g(x) — 0. But g(x) was the second derivative of log <p(x), hence log <p(x) 
is a linear function. Furthermore, log <p(.x) is periodic, which means that log <p(x) 
must be a constant. Therefore <p(x) is also a constant. We already know one value 
of ip(.v), namely, <p(0) = n. This implies that <p(x) — tt for all x. 

Recalling the definition of <p(x) in Eq. (4.1), we get 



n*)r (1 -x) = 



TT 

sin ti;.v ’ 



(4.5) 
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which is often called Euler's functional equation. If we set x = £ in Eq. (4.5), 
we get a new proof of Eq. (2.14). The exact value of the constant in Legendre’s 
relation was never used; therefore, this proof of Eq. (2.14) is independent of 
what was done in Chapter 3. 

With the help of Eq. (2.2) we can write Eq. (4.5) in the form 

n 

Sin ttx = r -= 7 r . 

- xr(x) i\- x) 

Now we substitute the expressions for J r ’(x) and T(— a) given by Weierstrass’ 
product formula, Eq. (2.8). This gives the following representation of sin 7r.xr as 
an infinite product: 

oo 2 

sin -nx = ttx 1 1 . (4.6) 

i=i v 1 

If we had assumed this product development of sin nx to start with (there 
are other proofs for the expression), we could have derived Eq. (4.5) by direct 
calculation. The approach we took is preferable, inasmuch as it gives us 
Eq. (4.6) at the same time. For the importance of Eq. (4.6) in analysis, we refer 
the reader to books on function theory. 

Let us formulate the main result of this section as a theorem: 

Theorem 4.1 

Every positive periodic function that has a continuous second deri- 
vative and satisfies the functional equation (4.3) is a constant. 
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Applications to Definite Integrals 



There is extensive literature that deals with more-or-less trivial manipula- 
tions of the two Euler integrals. We can only mention a few of the important 
results here. 

If we set e~‘ = r in Eq. (2.1), and then write t instead of t again, we get 

A*) = (I 08 1") 11 dU (5,l) 

Similarly, the substitution t x = r turns Eq. (2.1) into 

If we replace x by 1 /jc, then 

r(l + — ) = C e~'*dt. (5.2) 

\ X J J 0 

The case x = 2 in Eq. (5.2) is of special interest: 

C e^dt = ty;. (5.3) 

■' o 

The number of similar manipulations can be extended indefinitely. We 
will mention one other example of importance in analytic number theory. If 
a > 0, the substitution t — ar leads to the integral 

— ( f e- at t*~ l dt. (5.4) 

a* J o 

Now we turn our attention to the first Euler integral, Eq. (2.13). The 
substitutions t = t/(t + 1) and sin 2 <p , respectively, give the expressions 



and 



f x . r (x)r(y) 

Jo (I +t) x+u A* + y ) ’ 



(5.5) 



| (sin ipY x ~ l (cos <p) 2!/ ' 1 dtp = ^ 
' o 



r(x) r(y) 
A* + 3-) ' 



(5.6) 



28 



APPLICATIONS TO DEFINITE INTEGRALS 



29 



If we set y = (1 — x) in Eq. (2.13), Eq. (4.5) gives us the following inter- 
esting special cases: 



TT 


0 < x < 1, 


(5.7) 


sin ttx ’ 


TT 


0 < x < 1, 


(5.8) 


sin ttx ’ 


n 

sin ttx * 


0 < x < 1. 


(5.9) 



I £*-‘(1 - t)~ x dt 
o 

f* #»-» , 

dt 

J o 1 + < 



r” 1 * tt 

2 J (tg tpf x ~ X 4<p = -r— 

j n sin i 



If x andy are both rational numbers, Eq. (2.13) is the intergal of an algebraic 
function. Suppose we set x = mjn and y — \ in Eq. (2.13), and make the 
substitution t = t". We get 

JoVl-t” »r(mjn + 1/2) 

For m = 1 and n = 4 or n = 3, we get the following numerical results, with 
the help of Eqs. (3.11) and (4.5): 



dt 

Kvi~ 



(AjW 

? V327T 



f 1 dt (T(l/3)) 2 
Wl-f^ Vi V I6rr ’ 



which shows a connection between these particular numbers and the elliptical 
integrals. 

An integral representation for the error /u(.v) in Stirling’s formula can also 
be found. Since we have 

£7TT 2 '-<' + iM ,+ T)- 1 ' 

the series in Eq. (3.9) can be written in the form 

p ( *)=if 1 

’ “Vo* + ” + X 

If we define the following noncontinuous function: 

1% — t, for 0 < t < 1 

H(t) = 0, for t = 0 

\ otherwise periodic of period 1 
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it follows that 



„(*) = £ f A = £ C" J&L dt = lim f dt . 

+ n + x £ 0 '« * + x ’'-*J 0 t+x 



Because the integrand is an oscillating function that approaches zero as t 
approaches infinity, the integral 

(5.i3) 

exists. 

Equation (5.13) is the first step in deriving the so-called Stirling series, 
a refinement of Stirling’s formula. We now introduce the functions 



*U0 = 2(- 



cos 2int 

’ h w in 



«„.,(<) -2(-l)-|^V. (S.14) 

We know that — //,(/) = H(t), because — //,(/) is the Fourier series for //(f). 
For w 5= 2, the series H„(t) is absolutely and uniformly convergent for all x; 
for n — 1, the series converges uniformly in every closed interval that contains 
no integer. Therefore Eq. (5.14) implies that 

«;«(0 = UM (5.15) 

When n 2, this holds for all x ; when n = 1, for all nonintegral x. Because 
//„ + i(f) is always continuous, Eq. (5.15) implies that 

C H n (t) dt = //,,„(f) - H M {0). (5.16) 

J o 

The functions //„(<) are periodic of period 1, so it suffices to study them on the 
interval 0 ^ t < 1. Because H x (t) is a polynomial on this interval, H„(t ) must 
also be a polynomial there. We maintain that the coefficients of these polynomials 
are all rational numbers. This can be shown by induction. Our assertion is true 
when n = 1 ; we assume that it also holds for H n (t). This implies that 

f Unit) dt , 

J 0 

is a rational polynomial; thus it suffices to prove that //„ + ,(0) is a rational num- 
ber. When n + 1 is odd, Eq. (5.14) implies H n+1 (0) = 0. When n + 1 is even, 
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we can set H n+l (t) = 95(f) + H n+ j(0), where 95(f) is a rational polynomial. If we 
integrate this expression, then using Eq. (5.16) we get 

H n+i (t) - H n+i { 0) = // B+2 (f) = f‘ <p(t) dt + H n+l ( 0) /, 

J 0 

because n + 2 is odd. But Eq. (5.14) implies that // n+2 ( 1) vanishes; thus we 
have the expression 

0 = f l dt 4- H„ +l (0), (5.17) 

J 0 

which enables us to calculate H„ +i (0) and show that it is rational. 

We are now in a position to write Eq. (5.13) in the desired form. Repeated 
integration by parts gives 

’ x x 2 x 3 

, H„(0)(«-2)! r H n (t){n - 1)1 ^ tz \q\ 



rt B (0)(w -2)! _ r~ 

x"- 1 J 0 



(t + x)" 



An easy calculation shows that the sum of the last two terms in Eq. (5.18) is 
equal to 



f« ~ - l)\ 

Jo (' + *)" 



If rt is even, //„( 0) — // B (f) has the same sign (plus or minus) as H n ( 0) for all /; 
thus the integral above also has this sign. Furthermore, Eq. (5.14) implies that 
the numbers 7/ 2n+1 (0) all vanish, and that the signs of the numbers // 2 „( 0) 
alternate between plus and minus. This shows that the partial sums 

tf 2 (0)0! , tf 4 ( 0)2! , , H ln (0)(2n — 2)! 

T M I I* 



are alternately larger and smaller than /u(x). In other words, for every n there 
exists a number 6, 0 < 6 < 1 , so that 



f4*) = 



// 2 (0) 0! , H<( 0) 2! 



+ ••• + 



H^ 2 (0)(2n-4)\ 



+ — 2)!^ O<0<1. (5.19) 



We can not take the limit in Eq. (5.19) because the series diverges. But so 
long as n does not get too large, it gives us a very useful approximation. When 
n = 8, for example, we get the approximation 

/ \ _ 1 1 1 1 A (< 1Q\ 

_ 12x 360x 3 1260x6 1680x 7 - V ' 
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We now have a method for computing r(x). In the interval 4 < x < 5, 
for instance, we can compute /a(x) to within six decimal places by using Eq. (5.20). 

The first formula in Eq. (3.9) gives the value of log i T(.v) with the same accuracy. 

Using the functional equation (2.2), we finally get the value of F(x) in the interval 
from I to 2. This enables us to compute r(x ) for arbitrary x. When x is very 
large, Eq. (3.9) can be used directly. 

Anyone familiar with the so-called Bernoulli numbers will recognize the j 1 

connection with our H n ( 0). The reader is left to pursue this topic on his own. 
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Determining Y(x) by Functional liquations 



We have aquainted ourselves with three different functional equations for 
the gamma function: the functional equation (2.2), the multiplication formula, 
and Euler’s formula. To what extent is the gamma function determined by one, 
or a combination, of these equations ? 

Suppose for the time being that /(x) is an arbitrary function and that <p(x) 
denotes the quotient /(x)/.T(x). If /(x) satisfies a functional equation of the type 
in Eqs. (2.2), (3.10) or (4.5), then <p{x) will clearly satisfy the corresponding 
equation among the following: 

t <p( x + 1) = <p(*)> (6.1) 

<«> 

-*) = 1. (6.3) 

If /(x) satisfies Legendre’s functional equation (3.11), for instance, then 

9 (y) 9 (~y~) = <?>(*)• (6-4) 

For the sake of simplicity, we will assume from now on that f[x) satisfies 
the functional equation (2.2); consequently, <p(x) is periodic of period ], that 
is, <p(x) satisfies Eq. (6.1). We will also assume that /(x), along with <p(x), is 
continuous for all x. As a result of Eq. (6.1), the continuity of <p(x) for positive x 
implies continuity at zero and at negative integers, provided <p(x) is defined for 
these values in the right way. If Eqs. (6.2) or (6.3) also hold, they are valid for all 
x because of continuity. 

1 If we assume further that f(x) is always positive when .v is positive, then the 

logarithm of <p(x) is continuous. If we denote log <p(x) by g(x), the corresponding 
i functional equations (in addition to g(x + 1) = ^(x)) are 





> + ■ 


••+*( p ) =i{x) ' 


(6.5) 


g{x) = - 


«(1 


“*) = -£(-*). 


(6.6) 


*(t) 


+ g 




(6.7) 
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If f(x) has a continuous second derivative, then so does <p(x). We will 
assume this to be the case, and also that <p(x ) satisfies Eq. (6.4). But then <p(x) 
must be a constant. This follows from the theorem proven at the end of 
Chapter 4. Because of Eq. (6.4), the value of this constant must be 1. In other 
words, we have f(x) = C(.v). 

Theorem 6.1 

The gamma function is the only solution of both Eq. (2.2) and 
Eq. (3.11) that is positive for positive x and possesses a continuous second 
derivative. 



This theorem demonstrates the significance of the Legendre functional 
equation. Our next step is to find out whether the assumption of a continuous 
second derivative can be weakened. As a matter of fact, we will be able to show 
that a continuous first derivative is sufficient. 

Suppose we start with a preliminary observation. Equation (6.5) represents 
an infinite number of functional equations, one for each value of p. But these 
equations are not independent of each other. Assume, for instance, that Eq. (6.5) 
holds for the integers p l and p 2 . If we consider it for the integer p l , but with 
the argument (.* + k)jp 2 , we get 






Now we take the sum over k from zero to p 2 — 1 . On the right side we get g(x), 
since Eq. (6.5) holds for p 2 . But k + ip 2 runs over all integers from zero to 
PxPz ~ This yields the equation 



p,p a -i 

2 



2-0 



Plp2 ' 



therefore Eq. (6.5) also holds for the product p x p 2 . 

With this in mind, let us assume that Eq. (6.7) holds. Then Eq. (6.5) is 
valid for all integers of the form 2", and hence for arbitrarily large values of p. 
More generally, if Eq. (6.5) holds for an integer p, it also holds for the powers 
of that integer, and hence for certain arbitrarily large integers. 

Now we take the derivative of Eq. (6.5) 



P 



WtWFtV™ + *'( 



(X + p — 1 



)) = .?'(*)> 



( 6 . 8 ) 



and partition the positive x axis into intervals of length 1 ip, beginning with the 
origin. The arguments in the left side of Eq. (6.8) fall into p distinct consecutive 



intervals. Since g'(x) is periodic of period 1, we have a value of the function 
from each of the intervals between zero and I. But Eq. (6.8) holds for arbitrarily 
large values of p. As p approaches infinity, the left side of Eq. (6.8) converges 
to the integral 

fV(*)<fc«*(l)-j(0) = 0. 

J o 

This means that /(*) = 0 for all *; consequently g(x) is a constant. The value 
of this constant is zero, as can be seen from any of the equations in (6.5) which 
happen to hold. This proves the following theorem: 

Theorem 6.2 

The gamma function is the only continuously differentiable function 
that is positive for positive values of x, and that satisfies both Eq. (2.2) and 
Eq. (3.10) for some value of p. 

Now we might be tempted to conjecture that continuity alone suffices. 
This is not true at all. The function 

g( x ) = % 4 sin (2"™), (6.9) 

is continuous, because the series converges uniformly. It is also periodic, and it 
is easy to see that it satisfies Eq. (6.7). But it is not identically equal to zero. 
It is not even a constant: 



g( x ) = 0, g(l) = \ . 

If we assume mere continuity, what other properties must also be assumed 
to make Theorem 6.2 valid ? Equation (6.6) is not sufficient; the function 
in Eq. (6.9) also satisfies Eq. (6.6). A finite number of Eq. (6.5) is not sufficient 
either. Similar counter examples can always be constructed. But what happens 
if we assume that g(x) satisfies Eq. (6.5) for all integers pl As a conclusion to 
our study of the gamma function it will be shown that this property is sufficient. 

In order to do this we will make use of some facts about Fourier series. 
Let /(.v) be an integrable function of period 1 . Setting (for any a) 

c„ = ( f(x) e~ 2r,,n dx — [ f(x) e~ ini “ x dx, 

•' o •' a 

we associate with f(x) the series 




r «2 nivz 

> 
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regardless of whether or not this series converges. We denote this association by 

/(*) ~ X 

V= — CC 

We see immediately that the Fourier series of a sum of two functions is the sum 
of the two Fourier series. 

The Fourier series of f(x + a) has the coefficients 

/•I 

d y — f(x a) e _2m ‘' x dx = I f(x) ^-fc**^*-*) dx 
J o * 

fOt+l 

= e 4n<v * f(x) e- 2ni « dx; 

J a 

hence d v = e - n "' x c , . We see that the Fourier series of f(x + a) is obtained 
from that of f(x) by making the substitution x — »■ x + a. 

Let k ^ 1 be an integer. The Fourier coefficient d p of the function f(kx) 
is given by 

d p = f '/(**) dx = 4- (*/(•*) e~ 2ni ‘' x/k dx 

J o k J o 

= 4- X r f(x) e-‘ tnil ' xlk dx = 4- X fV(*) dx 

* S J tn-l * Tlx U 

= -j- (X f f(x) dx. 

k ,J o 

The sum in the last expression above is a geometric series with the 
ratio e~ 2nir,k . If v = fxk is divisible by k, each term of the sum is 1, and we 
obtain 

d„ k = f f(x) e -2 * 1 "* dx = c u . 

J A 



If v is not divisible by k, e~ tni,lk # 1 , and the formula for the sum of the geomet- 
ric series shows that d p = 0. Therefore 



+ <r, 

f(kx) ■ — ' 2) c t ,e >J,i,kx . 

fi =— ce 

The Fourier series for f(kx) is obtained from that of f(x) by merely substituting 
kx for x. 

If f(x) is continuous and the Fourier series converges for a particular value 
x„ , the value of the Fourier series for x = x 0 is the function value f(x 0 ). The 
reader will find a proof of this fact in most books on Fourier series. As a matter 
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of fact a consequence of Fejer’s theorem is that the Fourier series of a continuous 
function f(x) is always summable (even if not convergent) to f(x). 

We are now in a position to proceed with the proof at hand. Suppose our 
function g(x) satisfies Eq. (6.5) for all integers p. Replacing x by px, we obtain 
from Eq. (6.5) 

g( x ) + g (* + j) + ••• + g (* + “-) = g(P*)- (6.10) 



If we let 



g(x) ~ c p e tniyx . 



HI + 

g (x -| j ~ 2) C,e 2n "' z e inivml >’, 



g(p x ) ~ X c ' elt 



Substituting this into Eq. (6.10), we see that 



P-1 +x +ac 0-1 

X X c p e l7 ' i ' ,x e v ' i '' m l* = 2) c,e ini ' ,x e 2nivmlr '} 






is the Fourier series of the left side. 

Just as before, we obtain for 

v-\ 

^2 ntvm/p 

m=U 

the value p if v is divisible by p and zero if v is not divisible by p . The Fourier 
series of g(px) is, consequently, 

+•*> 

X fa* 

Comparing this result with Eq. (6.11), we find that 

C n = P C HV ■ 

In particular for the cases p = 1, — 1, 0 



c = £l • c — c — o 

l r> * > L -v ~ t, > c 0 — u. 
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Since we have assumed that Eq. (6.5) holds for all integers p ^ 1, we get 



CO 

g(x) ~ ^ e- m ' x -f- e -' lm%z ^ . 



But g(x) is a real-valued function, and hence 



CO ^ 

g(x) • — ^ sin 2vnx -f- — cos 2vnic^ . 



The terms 



^ sin 2m cos 2vnx 



are the Fourier series for the functions rr//(x) and — log (2 sin nx), respectively, 
where H(x) denotes the function introduced in Chapter 5. Therefore, at every 
point of continuity, 



g{x) = a-rrH(x) — b log (2 sin nx). 

The function H(x ) is bounded, the function — log (2 sin nx) is not. This 
implies that 6=0 for continuous g(x). But H(x) is not continuous; hence 
<t=0. This proves the following theorem: 



Theorem 6.3 

The gamma function is the only continuous function that is positive 
for positive x, and that satisfies Eq. (2.2) and Eqs. (3.10) for all values 
of p. 
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